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Abstract 



As recently found by Youm |hep-th/020126£ | , the entropy of the 



universe will no longer be expressible in the conventional Cardy-Verlinde 
form if one relaxes the radiation dominance state equation and instead 
assumes a more general equation of the form p = (7 — l)p, with 7 a 
constant. We show that Youm's generalized entropy formula remains 
valid when the cosmic fluid is no longer ideal, but is allowed to possess 
a constant bulk viscosity £. We supply our analysis with some numeri- 
cal estimates, thus calculating the scale factor a(t) for a k = universe, 
and also calculate via a perturbative expansion in £ the magnitude of 
the viscosity-induced correction to the scale factor if the universe is 
radiation dominated. 
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1 Introduction 



The suggestion of Verlinde |1J that there may exist a holographic bound 
on the subextensive entropy associated with the Casimir energy, makes it 
natural to ask: is this formal merging between the holographic principle, the 
Cardy entropy formula from conformal field theory ||, and the Friedmann 
equations from cosmology, merely a formal coincidence or does it reflect a 
deep physical property? As one might expect, the Verlinde proposal has 
inspired a considerable number of investigations, wherein various facets of 
the Cardy- Verlinde entropy formula, in particular, have been elucidated. An 
extensive list of references to the recent literature can be found in the paper 
of Youm || . Let us only mention here the large interest that has arisen in 
connection with the entropy and the energy as following from quantum and 
thermal fluctuations in conformal field theories |4j]. 

Verlinde assumed a spatially closed and radiation dominated universe. 
It becomes natural to inquire to what degree these basic assumptions can 
be generalized. For instance, Wang et al. || considered universes having 
a cosmo logical constant, and Nojiri et al. [[J considered quantum bounds 
for the Cardy- Verlinde formula. Recently, we examined the generalization 
of the Cardy- Verlinde formula allowing for a constant bulk viscosity ( in the 
radiation dominated universe |7|]. The conclusion of our analysis was that 
such a generalization is mathematically possible, although one should bear in 
mind that there occurs a possible interpretional difficulty: it is not evident 
from a physical point of view that the interpretiation of the subextensive 
part of the energy should be interpreted as a Casimir energy. 

The mentioned paper of Youm || is interesting, since it shows that the 
entropy of the universe can no longer be expressed in the conventional Cardy- 
Verlinde form if one relaxes the radiation dominance state equation and in- 
stead assumes a more general equation 

p=( 7 -l)p, (1) 

with 7 a constant. A more general form of the entropy formula results. And 
this brings us to the main theme of the present note, namely to explore to 
which extent the presence of a bulk viscosity in the cosmic fluid, together 
with the state equation (1), influences the entropy formula. In fact, it will 
turn out that the modified entropy formula found by Youm still persists, 
even when a constant bulk viscosity is allowed for. The central formula is 
Eq. (15) below. 
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We also supply the general analysis with some numerical estimates on 
viscous cosmology. Thus Section 3.1 calculates the scale factor a(t) for a 
viscous universe, assuming for simplicity that k = 0, if the viscosity is taken 
to set in suddenly at some instant t = t in , which we choose to lie in the 
plasma era. Section 3.2 exploits the smallness of ( to calculate the viscosity- 
induced correction to a(t) via a perturbative expansion in ( to the lowest 
order. The correction becomes very small, if ( is calculated in the plasma 
era from conventional kinetic theory. However, at earlier instants in the 
history of the universe the value of ( might have been much higher, thus 
strengthening the effect. 



2 Inclusion of viscosity 

We take the number n of spatial dimensions to be equal to 3, assume a FRW 
metric with k=+l, 

ds * = _ dt 2 + a 2 {t) ^_df_ + r 2 dQ 2^ ^ (2) 

set the cosmological constant A equal to zero, and assume that the cosmic 
fluid possesses a constant bulk viscosity (. The energy-momentum tensor 
can then be written as 



= P U^U U + (p - (0) V> ( 3 ) 

where in comoving coordinates £/ M = (1, 0), 9 = U^-^ = 3d/a, h^ v = g^ v + 
UpJJy being the projection tensor. The effective pressure is p = p — (6. From 
Einstein's equations — \Rg ilv = SirGT^ we derive the first Friedmann 
equation (the "initial value equation") 

H*= 8 -^p-\ (4) 



3 '" a 2 

(with H = a/ a), which is seen to be independent of the viscosity. The 
second Friedmann equation (the "dynamic equation"), when combined with 
Eq. (4), yields 

H = -4nG(p + p) + \, (5) 
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in which the effect from viscosity is explicit. The conservation equation for 
energy, T 0u . u = 0, yields 

P+(p + p)6 = (6 2 . (6) 
A direct calculation using the Friedmann equations leads to the equation 

^{pa i ) = (p-3p)a 3 . (7) 

So far, the equation of state for the cosmic fluid has not been used. Let us 
hereafter assume the specific form (1). Then, using Eq. (7) we derive 

f t (pa^)=(9W. (8) 

This equation is useful for our purpose. We wish to combine it with the 
conservation law for entropy: if n is the particle number density and a the 
entropy per particle, the four-divergence of the entropy current vector, = 
naW 1 , is 

S% = ^ (9) 
which, in view of the equation (nC/ M ) ;(U = 0, yields 

no = Ifl 2 (10) 

in the comoving coordinate system. 

We can now carry out the same kind of reasoning as in Ref . : Since ( is 
small, we can use for a = a(t) the same expression as for a nonviscous closed 
universe, namely a(t) = a* sin 77, where rj is the conformal time, the constant 
a* is defined as 

a* = ^(8nG/3)p in ai, (11) 

and the subscript "in" designates the initial instant of the onset of viscosity. 
(Strictly speaking, the approximate expression (11) assumes that the universe 
is radiation dominated.) Imagine that Eqs (8) and (10) are integrated with 
respect to time. Since (~ 1 pa 3j and C^na can be drawn as functions of t, it 
follows that pa 37 can be considered as a function of na. Then, since the total 
energy is E ~ pa 3 and the total entropy is S ~ naa 3 , it follows that Ea 3 ^' 1 ^ 
is independent of the volume V and is a function of S only. This generalizes 
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the pure entropy dependence of the product Ea, found by Verlinde |T[ in the 
case of a nonviscous radiation dominated universe. And it is noteworthy that 
the derived property of Ea 3 ^ -1 ^ formally agrees exactly with the property 
found by Youm || when £ = 0. 

Let us carry out the analysis a bit further, and write the total energy E 
as a sum of an extensive part E E and a subextensive part Eq'- 

E(S, V) = E E (S, V) + l -E c {S, V). (12) 

Under a scale transformation S — > XS and V — > AV with constant A, .E^ 
scales linearly with A. But the term i?c* scales with a power of A that is less 
than one: as Eq is the volume integral over a local energy density expressed 
in the metric and its derivatives, each of which scales as A -1 / 3 , and as the 
derivatives occur in pairs, the power in A has to be 1-2/3= 1/3. Thus we 
have 

E E (XS,XV) = XE E (S,V), E C (XS,XV) = A 1/3 E C (S, V), (13) 
which implies 

where C\,C<i are constants. In CFTs, their product is known: \fC\C2 = 
n = 3 [[J; this being a consequence of the AdS-CFT correspondence. Using 
Eqs. (12) and (14) we obtain 



S 



3 



E(j[2,E — E, 



c 



37-1 

(15) 



This is the generalized Cardy- Verlinde formula, in agreement with Eq. (20) 
in Youm's paper, reducing to the standard formula (with square root) in the 
case of a radiation dominated universe. In conclusion, we have extended the 
basis of Eq. (15) so as to include the presence of a constant bulk viscosity in 
the cosmic fluid. 



3 Numerical estimates 



5 



3.1 Spatially flat universe 

It is of physical interest to supplement the above considerations with some 
simple numerical estimates, showing, in particular, the order of magnitude 
of the viscosity terms. Let us assume, as mentioned above, that the effect of 
the bulk viscosity effectively sets in at some initial instant t = ti n and that 
( is thereafter constant for the times that we consider. For definiteness we 
choose 

t in = 1000 s (16) 

after the big bang. The universe is then in the plasma (or radiation) era; 
it consists of ionized H and He in equilibrium with radiation. The particle 
density is n in = 10 19 cm" 3 , and the temperature is T in = 4 x 10 8 K. The 
advantage of considering this relatively late stage of the universe's history is 
that the magnitude of ( can be calculated using conventional kinetic theory. 
At t = t in one finds 

C = 7.0 x 1(T 3 gcm^s" 1 (17) 

(cf. || and further references therein). 

For our estimate purposes it appears reasonable to assume that the in- 
fluence from the spatial curvature is not very important. For simplicity let 
us put k = 0. This implies that there exists the following simple differential 
equation for the scalar expansion ||: 

e(t) + ^9 2 (t)-l27rG(9(t) = Q, (18) 
which can be solved to give the expression for the scale factor: 

2/(3 7 ) 



a(t) = a- y 

where 



(19) 



t c = (UnGC)- 1 . (20) 
. The corresponding expression when viscosity is absent, is 

2/(37) 



a(t, C = 0) = a i: 



1 + ]pd- m (t - t h 



(21) 



The ratio between the expressions (19) and (21) reduces to unity in the limit 
when (t — t- m )/t c 1. This is we would expect. The influence from viscosity 
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generally turns up only in the factor t c , and the effect becomes strengthened 
when t c becomes smaller, i.e., when ( becomes larger. The numerical values 
given above for the instant t in = 1 000 s correspond to 

6 hl = 1 .5 x 10" 3 s _1 , t c = 5.1 x 10 28 s. (22) 



3.2 Perturbative expansion for a radiation dominated 
closed universe 

The smallness of ( makes it natural, as an alternative to the approach of the 
previous subsection, to make a Stokes expansion in (. For simplicity we now 
assume that 7 = 4/3, i.e., that the universe is radiation dominated. We put 
k — 1. Let subscript zero refer to quantities in the nonviscous case. We write 
the first order expansions as 

a = a (l + C«i), P = Po(l + CPi), (23) 

where the functions a x and pi are of zeroth order in (; they are regarded as 
functions of t or alternatively as functions of the conformal time r]. Corre- 
spondingly, the scalar expansion 6 = 3a/a is written as 

e = 6 (l + (6 1 ). (24) 

The Friedmann equation (4), and Eq. (8), now yield to the zeroth order 

61 = 24nGp - \, (25) 



a 



Pinaf n = PoOq, (26) 



and to the first order 



9<2i 

6 2 6 l = l27tGp oPl + ^, (27) 

A„af n (pi+4a 1 ) = 6 2 a 4 . (28) 
The solutions of Eqs. (25) and (26) are 

3 1 1 



a = a* sm 77, p = — — — 



SnG al sin 4 ?]' 
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#0 



3 COST} 

a* sin 2 r) ' 



where a* is defined by Eq. (11). From Eq. (28) we now have 
dpi 4da\ 



+ 



24nGa* sin rj cos rj, 



(29) 



(30) 



drj drj 

which can be integrated from rj = i] in onwards: 

Pi + 4ai = 87rGa*(cos 3 rj in — cos 3 77). 
We have here assumed that a\ = p\ = at the initial instant t = ti n . With 

(32) 



(31) 



a dai + 
Oi = —— tan 77 

drj 



(cf. Eq. (23)), we obtain from Eq. (27) a first order differential equation for 
the scale factor correction: 



sin 2t7 — — + 2(1 + cos 2 rj)ai + 87rGa* cos 3 ri = 8nGa* cos 3 % 
drj 



(33) 



This equation admits the integrating factor 

2(1 + cos 2 rj) 



Kv) = ex P 



sin 2rj 



drj 



V2 



sin 2 rj 

cos rj ' 



(34) 



whereby we find after some calculation the following solution, again observ- 
ing the initial conditions at t — t in : 



47rGa* 
sin 2 rj 



cos" // hl ( - cos2r/ — cos 2 ?7 in — - cos2r] in j cos // 



(35) 



Once the scale factor correction a% is known, the corresponding density 
correction p\ follows from Eq. (31). 

From Eq. (35) it is apparent that the relative correction £ai to the scale 
factor is of order 4irG(a^ or, in dimensional units, AirG^a^/c 3 . We here note 
that 4nG(/c 2 = l/(3t c ), thus about 6.5 x 1(T 30 s" 1 according to Eq. (20), 
which in dimensional form reads t~ l = \27rG(/c 2 . The quantity a*, according 
to Eq. (29), is the maximum value of the scale factor in a nonviscous k = 1 
theory. Let us put a* equal to the commonly accepted value of the radius 
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of the universe, i.e., a* = 10 28 cm. Then we obtain 47rG(a*/c 3 = 2 x 10~ 12 . 
The relative correction to the scale factor is thus in this case very small. 
Physically, this is due to the fact that we are considering an instant relatively 
late in the history of the universe. If the bulk viscosity ( were higher at earlier 
times, or, if there were an anisotropic stage present in the early universe at 
which the enormously higher shear viscosity would come into play ||, then 
the effect would be significantly enhanced. 
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